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Formulae for the energies of degenerate non-relativistic and ultra-relativistic Fermi gases play 
multiple roles in simple arguments related to the collapse of a stellar core to a neutron star. These 
formulae, deployed in conjunction with the virial theorem and a few other basic physical principles, 
provide surprisingly good estimates of the temperature, mass, and radius (and therefore also density 
and entropy) of the core at the onset of collapse; the final radius and composition of the cold compact 
remnant; and the total energy lost to neutrino emission during collapse. 



I. INTRODUCTION 
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Not long after 'supernovae' came to be distinguished 
from 'common novae' by virtue of their greatly higher in- 
trinsic luminosity^ and shortly after the discovery of the 
neutron in the early 1930s, Baade and Zwicky declared: 
"With all reserve we advance the view that supernovae 
represent the transitions from ordinary stars to neutron 
stars, which in their fin al s tages consist of extremely 
closely packed neutrons. For the 'core-collapse' va- 
rieties of supernovae (Types Ib/Ic/II) this basic picture 
prevails today with good theoretical and observational 
support)- 4 1 5 1 6 1 7 1 

The virial theorem and formulae for the energies of 
degenerate Fermi gases are especially important to this 
paper's simple arguments and calculations related to the 
collapse of the of a massive star's core to a neutron star. 
The virial theorem states that a self-gravitating body 
consisting of particles that are either non-relativistic or 
ultra-relativistic satisfies 



2.Ezv + Eji — —Eg, 
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where En and Er are the contributions to the body's in- 
ternal energy in non-relativistic and ultra-relativistic par- 
ticles respectively, and Eq is the gravitational energyP 
While hydrostatic equilibrium requires that pressure and 
density decrease with increasing distance from the center 
of a spherically symmetric star, an assumption of uniform 
density allows for simple and surprisingly good estimates. 
The gravitational energy of a uniform-density sphere of 
mass M and radius R is 



En = 
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where G is the gravitational constant. As will be seen 
below, degenerate gases dominate the internal energies 
of both the pre-collapse stellar core and the result- 
ing neutron star; hence formulae for the internal en- 
ergies (E F ) N and (E F ) R of non-relativistic and ultra- 
relativistic uniform-density degenerate Fermi gases con- 
fined to a spherical volume of radius R will play a promi- 
nent role in use ofEq. ([TJ. These energies are 



where Njy and Nr are the numbers of particles in non- 
relativistic and ultra-relativistic gases composed of a 
given particle type, and (cf)n and {zf)r are the respec- 
tive Fermi energies: 
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Here g is the spin degeneracy of the particle type, m is 
the particle mass (relevant to the non-relativistic case) , h 
is the reduced Planck constant h/2ir, and c is the speed of 
light (relevant to the ultra-relativistic case). In addition 
to their utility in representing the internal energies of 
self-gravitating bodies, Eqs. (|3|-((6| will also prove useful 
in a simple model of the nucleus used in arguments about 
the process of electron capture on nuclei that signals the 
onset of core collapse. 

After the existence of a degenerate iron core at the cen- 
ter of a massive star is argued for in Sec. [TTJ subsequent 
sections utilize Eqs. (|3])-((6| in various ways to determine 
the properties of this pre-collapse core (Sections III and 
IV ) and of the cold neutron star that is the final result 
of its collapse (Section [v}. The energy lost to neutrino 
emission — which allows collapse to proceed — is addressed 
in Section |VI| Section [VII| contains concluding remarks. 



II. A DEGENERATE IRON CORE 

Stellar collapse results from the evolution of a massive 
starPEl All stars spend most of their lives burning hydro- 
gen to helium in a central core, but in those with mass 
M > 8 M & (eight times the solar mass) the helium ash 
eventually ignites and burns to carbon, and then perhaps 
to oxygen and magnesium. For M > 11 M© additional 
core burning stages occur with increasing density and 
temperature, culminating with the production of iron- 
group nuclei; these lie near the top of the nuclear binding 
energy curve, and so no further energy can be extracted 
by nuclear fusion. 
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In very rough caricature, burning to higher-mass nu- 
clear species tends to proceed, at least at first, by suc- 
cessive captures of a particles (helium nuclei with pro- 
ton number Z = 2 and mass number A — 4, that is, 
|He). Two factors favor this. First, helium nuclei are 
unusually tightly bound, and this tends to make their 
presence in relatively high abundance energetically fa- 
vorable. Second, their modest proton number makes for 
a less formidable Coulomb barrier: it is much easier for 
an a particle to approach a higher-mass nucleus than it 
is for two higher-mass nuclei to approach each other. 

The dominant naturally-occuring isotope of iron is 
2@Fe, whose mass number A — 56 is indeed an integer 
multiple of 4. While it is not the most tightly bound nu- 
cleus overall,^ it is the most tightly bound nucleus near 
the valley of (3 stability,^ which tends towards greater 
neutron richness with higher A because of the increas- 
ingly high Coulomb energy of confined protons. 

The synthesis of ^Fc is not merely an unbroken se- 
quence of a captures. A variety of strong and electromag- 
netic interactions occur; these include fusion of higher- 
mass species when a particles are scarce, but the high 
temperatures required by the associated high Coulomb 
barriers also induce photodissociation that replenishes 
the a particle abundance, enabling the predominant role 
of species whose mass numbers are multiples of 4. More- 
over, the neutron excess manifest in 2 gFe requires the 
conversion of two protons to neutrons for each nucleus 
of this isotope, and so weak interactions must also occur 
along the way. 

Nevertheless, the assumption that the stellar core 
burns all the way to ItsFe, with a prominent role for a 
capture, allows the pre-collapse iron core temperature 
Tpo to be estimated. 

The first step in this estimate is to understand the re- 
lationship between temperature and the most favorable 
energy i^Gamow for a nuclear reaction to proceed, which 
is typically higher than kT (where k is Boltzmann's con- 
stant). This most favorable energy is the maximum of the 
so-called Gamow peak, which represents a compromise 
between the higher probability of Coulomb barrier pene- 
tration for higher energy and the declining number of par- 
ticles in the high-energy tail of the Maxwell-Boltzmann 
distribution! 4 ^ In particular -Ecamow is the maximum of 
the typically narrowly-peaked function eT b l E 1 e~ E / kT , 
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Here a » 1/137 is the fine structure constant (the two 
distinct uses of the symbol a — to denote the |Hc nucleus, 
or the fine structure constant — should be clear from the 
context); and 
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is the reduced mass of colliding nuclei of proton and mass 
numbers Z\,A\ and Z 2 ,A 2 , where tunc 2 = 931.5 MeV 
is the energy of the atomic mass unit. 

Normally the temperature is considered a 'given' and 
^Gamow is then determined from it; the strategy here 
will be to estimate the minimum value of Bo amow for 
which a capture to iron could readily occur, and then in- 
vert Eq. ^ to determine the corresponding temperature 
Tp . The uncertainty principle serves to characterize the 
point at which tunneling through the Coulomb barrier 
first becomes effective: 



Ax Ap c « Tic, 
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where the factor of c is included for calculational conve- 
nience (he = 197.33 MeV fm). Uncertainties in position 
Ax and momentum Ap are associated with the width 
and height of the barrier as illustrated in Fig. [I] The 
width of the barrier is 



Ax = r tp - r N I A 1 
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Here r tp is the classical 'turning point' separation at 
which the Coulomb energy equals i?Gamow) 
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and r^A^ 3 and r^A^ 10 are the radii of the colliding nu- 
clei; when they overlap, the deep energy well provided by 
the strong interaction dominates the potential governing 
the collision (r^ ~ 1.2 fm is the nucleon radius). 4 The 
momentum scale Ap associated with the height of the 
barrier above -Ecamow is given by 
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For a putative final a capture to ^Fe one has Z\ = 2, 

52. With these parameters, 
, = 8.7 MeV, 



Z 2 = 24, A l = 4, and A 2 
Eqs. yield r tp = 8.0 fm, Eq 

and kT Fc = 0.56 MeV. 

By this point the electrons in the iron core must be 
degenerate, thanks to neutrino emission. The thermal 
energy 



fcT Fe « 0.56 MeV 
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estimated above is comparable to the electron mass en- 
ergy m e c 2 = 0.511 MeV, so in the absence of an elec- 
tron chemical potential of similar magnitude positrons 
will be thermally produced. But these e + will anni- 
hilate with e~ into vv pairs that escape freely; hence 
we expect cooling by vv emission back down towards 



kT. 



0.5 MeV. In response to the resulting loss 



of pressure support the core contracts to higher density, 
with the conversion of gravitational potential energy in- 
creasing the temperature and bringing about further pair 
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FIG. 1: The potential governing a putative final a capture 
to jjjjjFe (heavy solid line) is taken to be Coulomb repulsion 
until the colliding nuclei begin to overlap, at which point the 
strong interaction provides a deep well. The minimum value 
of -EGamow for which tunneling is effective can be estimated 
with the uncertainty principle Aa; Ap w h, where Aa; is the 
effective width of the barrier (difference between the 'turning 
point' separation ftp and the 'overlap' separation), and Ap is 
related to the height of the potential barrier above Bcamow 
(jj, is the reduced mass of the colliding nuclei) . 



Use of Eq. ([15} and g = 2 in Eqs. ^ and Q for the 
core's relativistic electron energy En = (Ep) e in Eq. ([lj, 
together with E N w and Eq given by Eq. ([2]), results 
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The cancellation of radius is a clear indication that some- 
thing goes haywire when a star tries to support itself 
by relativistic degeneracy pressure alone. Equation (16 1 



nevertheless yields an estimate of the mass of this ex- 
treme configuration, the so-called Chandrasekhar mass: 



M Fo w 1.5 M G 



Y, 



Fc 



0.46 
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where Y e< p e sa 0.46 is the ratio Z/A of the proton and 
mass numbers of z^Fe. This may be taken as an es- 
timate of the mass of the core upon collapse, and of 
the compact remnant that results therefrom. (Taking 
inhomogeneous stellar structure into account yields 10 
M Fc «1.22Af Q (y e , Fc /0.46) 2 .) 



emission, and so on. Continued increase in density with 
temperature regulated to ~ 0.5 MeV/k would eventually 
result in electron degeneracy, which becomes the domi- 
nant source of pressure support. (Iron core central tem- 
peratures of ~ 0.5 — 1 MeV/fc are in fact seen in stel- 
lar evolution codes, though the neutrino pair emission 
that eventually gives rise to degeneracy actually begins 
as early as carbon burning!^ when e + e~ pairs begin to 
be produced from the high-energy tail of the photon dis- 
tribution.) 



III. MASS OF THE CORE 

Catastrophic collapse ensues when the iron core be- 
comes sufficiently massive that it can no longer be sup- 
ported against gravity by electron degeneracy pressure. 
As the mass and density of the iron core increase dur- 
ing the final burning stages, the velocity of the degen- 
erate electrons approaches the speed of light; after this, 
the electrons' contribution to pressure support cannot be 
further increased. 

An estimate of the core's maximum mass follows from 
consideration of a uniform core whose pressure support 
comes solely from a degenerate Fermi gas of relativistic 
electrons. The number of electrons in the core is 



N K = 
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where M now denotes the mass of the core rather than 
the mass of the entire star, and the electron fraction Y e 
is the ratio of the net number density of electrons n e - — 
n e + (or, by charge neutrality, the total number density 
of protons n p ) to the total baryon number density n. 



IV. DENSITY AND ENTROPY PER BARYON 
OF THE CORE 

The vigorous onset of electron capture — one proximate 
cause of the core's instability as it approaches the Chan- 
drasekhar mass — can be used to estimate the density of 
the iron core just before collapse. This new physical in- 
put is needed because, in a uniform-density model, the 
density follows from knowledge of the radius R, which 
cancelled out of the estimate of the core mass M in 
the previous section. Electron capture e~ + (Z, A) — > 
(Z — 1, A) -\-v e begins in earnest when the electron Fermi 
energy 



1/3 

V 4 m N ) R 



/9tt MY e \ 1/J he 



(18) 



rises above the energy threshold required to convert a 
proton to a neutron — when both of these baryons are 
bound in a nucleus. 

Hence calculation of the energy threshold for electron 
capture requires some model of the nucleus; fortunately, 
the concept of degeneracy once again facilitates a sim- 
ple estimate. Crudely considering a nucleus of proton 
number Z and mass number A as consisting of degener- 
ate Fermi gases of non-relativistic protons and neutrons 
confined to a volume of radius ta = r^A 1 ^ 3 , the typ- 
ical emitted neutrino energy e„ e is not (ei?) e , nor even 
(ep) e — A, where A = (m n — m p — m e )c 2 » 0.78 MeV is 
the difference between the neutron mass energy m n c 2 and 
the proton and electron mass energies m p c 2 and m e c 2 : in 
addition to A, the energy required to lift a bound pro- 
ton from its Fermi sea (within the nucleus) to the top of 
the neutron Fermi sea (within the nucleus) must first be 
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paid. Hence e„ e w (e_F) e - A - [(e_F) n;j4 ,z - (e_F) P;J 4,z], 
where (e F ) n -,A,z and (e F ) n] A,z are the Fermi energies of 
the neutrons and protons bound in the nucleus. Using 
Eq. |5]) for these Fermi energies results in 
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where it is assumed that all baryons are locked up in 
nuclei characterized by the same Z and A. 

Electron capture can proceed when the right-hand side 
of Eq. (20 1 is positive; using (e F ) e as given in Eq. (18) 
yields Rp c « 8.1 x 10 2 km at the onset of electron capture 
for Mp e = 1.5 Mq and Ve.Fe = 0.46, corresponding to 



p Fc w 1.3 X 10 12 kg/m 3 
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in the iron core at the onset of collapse. 

Knowledge of the density then enables the entropy per 
baryon of the core to be estimated: sf c ~ 0.91 k receives 
contributions from iron nuclei of mass number A = 56, 
relativistic degenerate electrons, and photons: 
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where vtla ~ m^A is the mass of a nucleus of mass 
number A, fcTp c ~ 0.56 MeV as found in Sec. |TTJ and 
(ei?) e ,Fc ~ 4.4 MeV is the electron Fermi energy obtained 
fromEq. (fill. 



(As was the case with Tp e , the values of pFe and SFe ob- 
tained here are in reasonable agreement with the central 
iron core values produced by stellar evolution codes.^) 



V. COMPOSITION AND RADIUS OF THE 
COMPACT REMNANT 

Having characterized the precollapse core and deter- 
mined that it is doomed to collapse — and supposing that 
the degeneracy pressure of nonrelativistic nucleons will 
halt collapse — what might we expect the final radius 
-Rcoid of the compact remnant to be, and why do we ex- 
pect it to consist of neutrons rather than a mixture of 
neutrons, protons, and electrons? 

Aware that weak interactions can interchange neutrons 
and protons, we provisionally take our star of "extremely 



closely packed [nucleons]" to consist of nonrelativistic de- 
generate Fermi gases of neutrons and protons and a rel- 
ativistic degenerate Fermi gas of electrons. The energies 
(Ep) n and (E F ) p of the nucleons are 
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according to Eq. ([5|. The virial theorem of Eq. (fij, 
with E N = (E F ) n + (E F ) p , and E R = (E F ) e and E G 
as these appear on the left- and right-hand sides of Eq. 
( 16 1 respectively, provides one equation describing the 



cold compact remnant. A second constraint is provided 
by chemical ('/?' or 'weak') equilibrium, which requires 



(29) 



where the chemical potentials fii include particle masses. 
That is, Hi = rrii + (e F )i at zero temperature, where (e F ) n 
and (e F ) p are given above and (e F ) e is given by Eq. ( 18 1. 
Taking M = Mp e = 1.5 Mq, simultaneous solution of 
Eqs. (fij) and (29) as specified here yields i? m id = 11 km 
and Ye, cold = 1.0 x 10~ 2 . Because Y e is the fraction of 
nucleons that are protons, the label 'neutron star' is well 
deserved indeed! This neutronization occurs (via elec- 
tron capture) because the increasing density associated 
with collapse makes the degeneracy energy of relativistic 
electrons increasingly costly. 

Neglecting the presence of protons and electrons — that 
is, taking Y e ss — the virial theorem becomes (E F ) n — 
—Eq/2, and this alone yields 



R, 



cold 



11 km 



1.5 M G 
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(The resulting baryon number density 



0.25 fin" 



is within a factor of two of the value (47rr^/3) 1 ~ 
0.14 fm~ 3 derived from from experimental measurements 
of the size of atomic nuclei, which also feature densely- 
packed nucleons. In some respects, a neutron star is not 
unlike an atomic nucleus of astronomical proportions!) 



VI. ENERGY LOSS TO NEUTRINOS 

The vigorous onset of electron capture not only signals 
the onset of collapse, as discussed in Sec. |IV[ the energy 
loss by neutrino emission incident to continued electron 
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capture is necessary for collapse to continue to comple- 
tion. This is another consequence of the virial theorem: 
from Eq. ([I]) we see that the total energy of the prec- 
ollapse core, whose internal energy is completely domi- 
nated by relativistic electrons, is (Ep) e + Eq w 0; while 
the total energy of the cold neutron star, dominated by 
nonrelativistic neutrons, is {Ep) n + Eq w Eq/2. Hence 
an amount of energy equal in magnitude to half the final 
neutron star gravitational energy, 

*-"*»"'G^)"(^)- w 

must be lost to neutrinos (and a bit of it to the supernova 
explosion) before the neutron star can reach its cold final 
state. 



VII. CONCLUSION 

The observed supernova phenomenon, which led Baadc 
and Zwicky to speculate on the formation of neutron 
stars, exhibits explosion energies of about 10 44 J — a tiny 
fraction of the 10 46 J of gravitational potential energy 
released in stellar collapse. Because the explosion is a 
comparatively minor detail in energetic terms, it is not 
surprising that sophisticated simulations are required to 
understand it! 11 * 12 ! 

Nevertheless, some features of the stellar core collapse 
itself can be understood very simply, using estimates 
that make use of the virial theorem and formulae for 
the energies of degenerate gases. (More detailed but 
less self-contained studies of the nuclear physics, weak 
interactions, and hydrodynamics related to the state of 
the pre-collapse core and the process of collapse, which 



retain an analy tic or se mi-analytic flavor, are available 
in the literature! 13 * 14 * 15 ^) The virial theorem, which re- 
lates internal energy to gravitational energy, is a concise 
summary of the hydrostatic equilibrium that determines 
stellar structure; an assumption of uniform density then 
provides simple approximate relations between the mass 
and radius of both the pre-collapse stellar core and the 
post-collapse neutron star. The pre-collapse core is sup- 
ported by electron degeneracy pressure, and when these 
electrons are relativistic an estimate of the core mass 
just before collapse is obtained. The radius (and there- 
fore density) of the pre-collapse core can be obtained by 
estimating conditions for the onset of electron capture, 
through which energy (and pressure support) are lost 
through neutrino emission. This calculation involves an- 
other use of the concept of degeneracy, namely to crudely 
model the nucleus as consisting of degenerate proton and 
neutron ideal gases. Taking the cold post-collapse neu- 
tron star to be supported by degenerate non-relativistic 
nucleons, the virial theorem and the condition of 'chem- 
ical' equilibrium provide an estimate of the neutron star 
radius and demonstrate why the name 'neutron star' is 
well-deserved. Finally, the difference in total energy of 
the pre-collapse core and the post-collapse neutron star 
implied by virial theorem constitutes an estimate of the 
energy loss to neutrino emission that allows collapse to 
occur. 



Acknowledgments 

This work was supported by Oak Ridge National Lab- 
oratory (ORNL), managed by UT-Battelle, LLC, for the 
DoE under contract DE-AC05-00OR22725. 



* cardallcy@ornl.gov 

1 W. Baade and F. Zwicky, "On Super-novae," Proc. Nat. 
Acad. Sci. 20, 254-259 (1934). 

2 W. Baade and F. Zwicky, "Supernovae and Cosmic Rays," 
Phys. Rev. 45, 138 (1934). 

3 W. Baade and F. Zwicky, "Cosmic Rays from Super- 
novae," Proc. Nat. Acad. Sci. 20, 259-263 (1934). 

4 D. Arnett, Supernovae and Nucleosynthesis (Princeton 
University Press, Princeton, 1996). 

5 S. E. Woosley, A. Heger, and T. A. Weaver, "The evolution 
and explosion of massive stars," Rev. Mod. Phys. 74, 1015- 
1071 (2002). 

6 A. V. Filippenko, "Optical Spectra of Supernovae," Annu. 
Rev. Astron. Astrophys. 35, 309-355 (1997). 

7 J. M. Lattimer and M. Prakash, "The Physics of Neutron 
Stars," Science 304, 536-542 (2004). 

8 G. W. Collins II, The Virial Theorem in Stellar As- 
trophysics, vol. 7 of Astronomy and Astrophysics Series 
(Pachart Publishing House, Tucson, 1978). 

9 M. P. Fewell, "The atomic nuclide with the highest mean 
binding energy," Am. J. Phys. 63, 653-658 (1995). 



S. Chandrasekhar, "On stars, their evolution and their sta- 
bility," Rev. Mod. Phys. 56, 137-147 (1984). 

11 A. Mezzacappa, "Ascertaining the Core Collapse Super- 
nova Mechanism: The State of the Art and the Road 
Ahead," Annu. Rev. Nucl. Part. Sci. 55, 467-515 (2005). 

12 S. Woosley and T. Janka, "The physics of core-collapse 
supernovae," Nature Physics 1, 147-154 (2005). 

13 H. A. Bethe, G. E. Brown, J. Applegate, and J. M. Lat- 
timer, "Equation of state in the gravitational collapse of 
stars," Nucl. Phys. A 324, 487-533 (1979). 

14 G. M. Fuller, W. A. Fowler, and M. J. Newman, "Stel- 
lar weak interaction rates for intermediate-mass nuclei. II. 
A=21 to A=60," Astrophys. J. 252, 715-740 (1982). 

15 G. M. Fuller, W. A. Fowler, and M. J. Newman, "Stellar 
weak interaction rates for intermediate-mass nuclei. IV. In- 
terpolation procedures for rapidly varying lepton capture 
rates using effective log(ft)-values," Astrophys. J. 293, 1— 
16 (1985). 

16 A. Yahil, "Self-similar stellar collapse," Astrophys. J. 265, 
1047-1055 (1983). 



